The Polyakov loop and the heat kernel expansion at finite temperature 
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The lower order terms of the heat kernel expansion at coincident points are computed in the 
context of finite temperature quantum field theory for flat space-time and in the presence of general 
gauge and scalar fields which may be non Abelian and non stationary. The computation is carried out 
in the imaginary time formalism and the result is fully consistent with invariance under topologically 
large and small gauge transformations. The Polyakov loop is shown to play a fundamental role. 
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The heat kernel e^"^^, where K is an Euclidean Klein- 
Gordon operator, was introduced by Schwinger in 
quantum field theory as a tool to regularize ultraviolet 
divergences in a gauge invariant way. The procedure of 
systematically isolating these divergences was developed 
by DeWitt through an asymptotic expansion in powers 
of the proper time r (2[ ||] . The heat kernel as well as its 
expansion has also been applied to study spectral den- 
sities and indices of Dirac operators {D) ^ in terms 
of effective Klein-Gordon operators to compute 

the function Q and anomalies of such operators , to 
properly define the effective action of chiral gauge theo- 
ries to the Casimir effect Q, etc. Exact calculations 
of the heat kernel at coincident points are available in 
particular manifolds 0, and the coefficients of the 
expansion have been computed to rather high orders in 
several setups, including curved spaces with and without 
boundary, and in presence of non Abelian gauge fields 
and non Abelian scalar fields, using different methods 
|, ^, In this work we deal with the 

extension of these results to finite temperature in flat 
space-time. 

As is well known, in the Matsubara formalism for quan- 
tum field theory at finite temperature, wave functions are 
required to be periodic for bosons and antiperiodic for 
fermions, with period /3 (the inverse temperature) and 
the external fields are also periodic. This introduces a 
compactified Euclidean time . The heat kernel can be 
regarded as the evolution operator (in imaginary proper 
time r) of a first quantized particle [Q, thus it admits 
a path integral representation, which at finite tempera- 
ture must include also multiply connected paths in the 
compactified space-time. The standard results quoted 
e.g. in 1^ for the heat kernel expansion refer to general 
smooth manifolds without boundary and thus, in prin- 
ciple, they would directly apply to finite temperature as 
well. However, the standard formulas are not sensitive to 
the compactification of time; they only sum paths which 
are topologically trivial (contractile to a point). Actually 
this correctly reproduces the analytical part in a strict ex- 
pansion in powers of the proper time, but neglects the ex- 
ponentially suppressed finite temperature contributions 



]l9[ . To see any temperature effect in the heat kernel 
expansion it is necessary to include also the paths that 
wind once or more around the compactified time. This 
is automatically achieved in the canonical formalism by 
replacing the continuous energy variable by the discrete 
Matsubara frequencies. More generally, as already ar- 
gued in (see also (2|]), in the presence of gauge fields, 
finite temperature effects are tied to the presence of the 
(untraced) Polyakov loop 



Vt{x) = Texp - 



(1) 



(T refers to temporal ordering, this is required to allow 
for the non Abelian case). The presence of ^{x) is to be 
expected on general grounds. Let us regard the heat ker- 
nel at finite temperature as the evolution operator of a 
first quantized particle. The space-time manifold defines 
a d-dimensional cylinder which plays of role of space for 
the particle and —it is the time. If the circulation of the 
"vector potential" around the cylinder is non van- 
ishing this signals the presence of an Aharanov-Bohm 
magnetic flux through its interior, which will affect the 
propagation. The particle will pick up a phase f2(a;) after 
looping once around the cylinder (in the general case a 
non Abelian phase |^). 

A computation of the thermal heat kernel coefficients 
was presented in [p3| where a common renormalization 
for all the coefficients, — /(T//3^)a„, was found. How- 
ever, because the temperature effects account for the fact 
that Euclidean time is compact, an explicit breaking of 
Lorentz invariance is expected and the previous factor- 
ization cannot hold in general (see e.g. The former 
computation refers to flat space-time (without bound- 
aries) and non Abelian scalar and gauge flelds. In 
a similar study is presented where curvature is allowed 
but fields are assumed to be static. These authors find 
that /(r//3^)a„ is augmented with new terms that break 
Lorentz symmetry, so-called anomalous terms also found 
in However both calculations are incomplete since 
they are in conflict with results for Abelian and station- 
ary configurations (in flat space) presented in pj| . There 
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it is shown that the Euclidean scalar potential Aq {x) can 
appear in the heat kernel (outside covariant derivatives) 
provided that it does so as a periodic variable with pe- 
riod 27r//3, to be consistent with (large) gauge invariance. 
This is a particular manifestation of the relevance of the 
Polyakov loop. 

We present here the first systematic calculation of the 
heat kernel expansion at finite temperature at coinci- 
dent points (in flat space without boundaries) for general 
backgrounds which may be non Abelian and non station- 
ary. The result is gauge covariant and depends locally 
on the covariant derivatives and the Polyakov loop. The 
Klein-Gordon operator to be considered is of the type 
K = M — Df^, where = 9^-|-A^(a;) is the gauge covari- 
ant derivative and M{x) is a Lorentz scalar. Both A^(a;) 
and M{x) are allowed to be matrices in some internal 
space (the treatment applies to Abelian and non Abelian 
theories) . One can form local and gauge covariant opera- 
tors by combining the basic operators M, F^^, — [I?^, D^] 
and their covariant derivatives, e.g. O = F^^Vi^J^^M , 
where V^^X stands for The local operators can 

be classified according to their mass dimension. For 
each dimension there is a finite number of them. At 
zero temperature, under suitable regularity conditions, 
one can asymptotically expand an expression of the form 
{x\f{K)\x) as a sum of local operators, i.e. 

(a;|/(i^)|a;) -^/,a(a:), (2) 

r 

Or denotes a linear basis including all possible local op- 
erators and fr are some numerical coefficients which de- 
pend on the function /. In particular, for the matrix 
elements of the heat kernel at coincident points one ob- 
tains the familiar Seeley-DeWitt asymptotic expansion 

oo 

(xle-^^lx) = (4^t)-'^/2 ^ a„(a;)r" , (3) 

d being the space-time dimension. The local operators 
Unix) themselves do not depend on cZ or r and have mass 
dimension 2n. Of course, once one has the expansion for 
the heat kernel, the use of integral representations |l[] al- 
lows to obtain the expansion for other functions f{K), 
such as the ^-function ^, the effective action, the 
resolvent or the effective current, without a new compu- 
tation. 

At finite temperature, because the Euclidean time 
coordinate is compactified, one can expect to expand 
{x\f{K)\x) in terms of operators which are local regard- 
ing the spatial coordinates but not in the time coordinate 
^7| . Nevertheless, the following asymptotic expansion 
still holds |§ 

{x\fiK)\x}=J2frm^))OAx). (4) 



That is, at finite temperature the numerical coefficients 
fr become functions of the Polyakov loop fl{x). This re- 
sult is remarkable on several accounts. First, although 
full space-time locality does not hold at finite temper- 
ature, still the result looks space-time local if written 
in terms of the field n{x) and the other local opera- 
tors. (Note that n{x) is local regarding gauge transfor- 
mations.) Second, full gauge covariance (the covariance 
under continuous local similarity transformations in in- 
ternal space) is manifest, since this property is shared by 
the Polyakov loop and the local operators. This is not 
always automatically true for other expansions; because 
the space-time is not simply connected in the imaginary 
time formalism, in addition to the usual topologically 
small gauge transformations there are also large transfor- 
mations (i.e., not reachable by composition of infinitesi- 
mal ones). Covariance under such large transformations 
fails to hold, for instance, in simple-minded perturbation 
theory ||^, |2||. 

In the non Abelian case r2(x) needs not commute with 
the other quantities. The fact that it appears to the left 
in (U) is a matter of choice (putting it to the right would 
yield a different set of functions fr). We remark that the 
operators Or are local and gauge covariant but not neces- 
sarily Lorentz invariant. That is, operators with different 
Lorentz indices may have different weights at finite tem- 
perature. Rotational invariance is of course preserved. 

As we have argued above, for many purposes it is suf- 
ficient to consider the expansion for the heat kernel. At 
finite temperature the result can still be written as 

(a;|e--^|a:) = (47rr)-''/2^aJ(x)r". (5) 

n 

The thermal heat kernel coefficients contain local op- 
erators, and depend on f2, r and /3, but not on d. 

At finite temperature the heat kernel expansion con- 
sists actually of three expansions. One in powers of M , 
another in powers of Vi, and finally another expansion 
in powers of Vq. The first two expansions are standard 
and no different to those done at zero temperature, but 
the expansion in Vq is tricky and has to be carried out 
with care. This is because a naive expansion in tempo- 
ral derivatives is in conflict with gauge invariance. This 
is easy to understand by noting that the derivative ex- 
pansion refers to external fields with very wide profiles. 
Whereas this is unproblematic in the spatial direction, 
the situation in the temporal direction is different since 
the Euclidean time is periodic and dilatations in the tem- 
poral direction break the periodicity condition. 

For the calculation we use the method of symbols 
which is based on the identity ^ 

{x\f{D,,M)\x)^^J2 J ■^;^{^\fiDf.+^P,,M)\0). 
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f{Dfj_, M) is any (ultraviolet finite) operator constructed 
out of and M, and |0) is the zero momentum state, 
(a;|0) = 1. In practice this means that 9^ in vanishes 
when it reaches the right end, and the brackets {x\ ■ \0) are 
often not explicit. An application of the symbols method 
formula in the temporal direction only, gives 



(x e 



Q = ipQ + Do 



{M~Q'^~D'i) 



(7) 



the replacement VI — f2. The (/3„'s are related to the 
special functions 6*3 and 9i and their derivatives ||3|] . 

The first thing to note is that, due to the summation 
over Po, the VJ„(r2)'s are one- valued functions of f2 and 
so gauge invariance is preserved. (The quantity logfi is 
covariant under small gauge transformations but many- 
valued under large ones.) Actually, if we use Poisson's 
summation formula for, say, the function tpo we obtain 



(bosonic) 



(12) 



For the spatial part we can directly use the known 
heat kernel expansion at zero temperature 12 on 
the effective (d — l)-dimensional Klein-Gordon operator 
(M — Q^) — The operators Q are then moved to 
the left generating commutators of the form [Q, ]. These 
commutators are just T>o = [Dq, ] since po is a c- number. 
The Q's not in commutators, placed at the left of the 
expressions so obtained, come in the form Q"e^'5 (cf. 
(O)). The final step is to use the operator identity pM 



(8) 



Because all external fields are periodic, e~^^° ~ 1 and 
g-/3-Do ig equivalent to Q{x). In this way the operator Q 
defined in (0) reduces to a function of x 



(9) 



An explicit computation gives, up to operators of mass 
dimension four [BTI, 
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In these formulas Ei = Foi is the electric field and we use 
the convenient notation M^^ = Vf^V^M, E^i = V^Ei, 
etc. The quantities (pn are functions of Vl{x) and the ratio 
tI0^, namely 



(47rr) 



1/2 



/3 



(11) 



Here, Q given in (||) and po runs over the Matsubara 
frequencies, 27rfc//3 [k integer) for bosons and 27r(fc -|- 
for fermions. Consequently, there are bosonic and 
fermionic versions for each function (^„, related through 



Second, we note that the expansion contains coefficients 
of half-integer order with local operators of odd dimen- 
sion. Due to the property (/'ri(^^) = (— l)"Vn(^~^) the 
symmetry under Do — > —Dq, which is manifest in the 
Klein-Gordon operator, is preserved. Odd dimensional 
terms vanish at zero temperature except in presence of 
boundaries One of the techniques to deal with 

boundaries is the method of images. The path integral 
method suggests that the presence of these terms at finite 
temperature can also be understood from the fact that 
the particle sees its own image from its periodic replicas. 

In the zero temperature limit the functions ipn become 
constant. For even n, they take the values (— ^)"^^(« — 
1)!! and vanish for odd n. It is easy to verify that the 
thermal coefficients reproduce the zero temperature ones 
in this limit and Lorentz invariance is restored. As illus- 
trated in (|l^) the low temperature or low t corrections 
are exponentially suppressed in a J. 

The thermal coefficients have been collected according 
to the mass dimension of the local operators. At zero 
temperature this is equivalent to an expansion in powers 
of T, but at finite temperature it corresponds to expand 
in powers of t keeping t//3^ and fl constant. Due to the 
exponential suppression, a strict expansion in r for fixed 
P would yield just the zero temperature result. (This 
accounts for the well known fact that ultraviolet diver- 
gences and anomalies are temperature independent p3| .) 

Also useful is the expansion of the trace of the heat 
kernel operator Tr(e~'^^) = J2n / d'^xtr{b'^). To mass 
dimension six we find 



bl = 



hi = -^oM 
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ipi{M„E,} 
1 



(13) 



6' 



12 



-^ip2E,ME, - (j^Vo + ^(^52 
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As is well known, the expression for 6^ is intrinsically not 
unique due to integration by parts and the trace cyclic 
property. At zero temperature these operations do not 
mix different orders, but at finite temperature this is no 
longer true since the commutators (in particular covari- 
ant derivatives) of Q,{x) generate dimensional terms H 



This explains why, e.g. 63^2 vanishes although Tr(a3^2) 
is not identically zero. Of course, the trace of the heat 
kernel itself is unique and the ambiguity amounts to a 
reorganization of the asymptotic expansion (no such am- 
biguity exists for the a^). The heat kernel is symmetric 
under transposition ABC ■ ■ ■ — > • ■ • CBA and this sym- 
metry is manifest in the (using the cyclic property). 
It can be noted that, exactly as in the zero temperature 
case [H, the coefficients aJ^ can be obtained from the &^ 
through functional derivation with respect to AI{x). 

In addition to the explicit formulas, one of the inter- 
esting results of this calculation is of qualitative type, 
namely, the operators appearing in the coefficients or 
&^ can be grouped according to the weights (fk- This in- 
dicates a definite pattern of Lorentz symmetry breaking 
at finite temperature (even in the absence of scalar po- 
tential). This pattern will manifest itself in applications 



sucH as those to renormalization group how at hnitc tem- 



perature |p5[, the general form ot the chiral perturbation 
theory effective Lagrangian at finite temperature, etc. 

The fundamental role played by the Polyakov loop is 
further illustrated by considering the trivially solvable 
case of a homogeneous relativistic quantum gas. In the 
absence of external fields (other than a mass to) , the low- 



tive I action and so the grand canonical potential 



exactly factorizcs out of the heat kernel) A simple 
integral transform to pass from the exponential to the 
logarithm function reproduces the standard result, with 
distribution function (e''('^'»-*') - l)-i - (e'^('^>'+^'^ - 1)-^ 
(for bosons) [Q. We want to emphasize the relation be- 
tween the chemical potential fi and the Polyakov loop 
il. The chemical potential couples as a constant addi- 
tive term in the scalar potential. Because it is constant, 
fj, does not contribute to the local operators, since there 



Aoi'i:) only appears Llii'ough lha covai'laiiL dci'lvatlvu Vq. 
If the Polyakov loop were abseuL m Llie fomiulas, fi would 
not appear at all in the partition function. An obviously 
inco rrect result. It is also noteworthv that the periodic 



dep endence of the heat kernel on log CI implies the well 



known fact that the partition function is periodic in 
with period 27ri (a consistency condition due to its cou- 
pling to the quantized charge operator). The Polyakov 
loop thus appears as a generalization of the factor e^^ 
to non constant and non Abelian gauge fields. In curved 
space-times, in addition to the Polyakov loop of the gauge 
connection , there is a Polyakov loop tied to the par- 
allel transport connection P^, with consequences in field 
theory in presence of gravitational fields. 
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